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Abstract 

We study a special case of the Brauer model in which every path of 
the model has weight q = 1. The model has been studied before as a 
solvable lattice model and can be viewed as a Lorentz lattice gas. The 
paths of the model are also called self-avoiding trails. We consider the 
model in a triangle with boundary conditions such that one of the trails 
must cross the triangle from a corner to the opposite side. Motivated 
by similarities between this model, SLE(6) and critical percolation, we 
investigate the distribution of the hull generated by this trail (the set of 
points on or surrounded by the trail) up to the hitting time of the side 
of the triangle opposite the starting point. Our Monte Carlo results are 
consistent with the hypothesis that for system size tending to infinity, 
the hull distribution is the same as that of a Brownian motion with 
perpendicular reflection on the boundary. 

1 Introduction to the model 

In this paper we present results from a Monte Carlo study of a special case 
of the Brauer model. This model has appeared in the literature in different 
guises, and received the name Brauer model only recently. Originahy, it 
was studied as a g-state solvable vertex model 120] and later as an 0{q) 
symmetric, solvable lattice model l2IS j . More recently it was observed 
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that the model could be seen as a model of paths on the lattice in which 
each closed path has a weight equal to q, where q can take non- integer 
values. In this language of paths each vertex can carry one of the following 
configurations of path segments, with the corresponding weights: 



Here the weights are chosen to be solutions of the Yang-Baxter (or star- 
triangle) equation [31 . Although the integrability condition in the Yang- 
Baxter equation does not restrict q, the actual solution by means of the 
Bethe Ansatz was only extended to the integers |17| . 

The name Brauer model was given to emphasise that the transfer matrix 
of the model is an element of the Brauer algebra [Hj- The model attracted 
particular attention in the limit g — 1 when it describes the probability 
distribution of lattice paths taken by a particle that is scattered by randomly 
placed, static scatterers [SJ^lElUn]- As such it can be viewed as a Lorentz 
lattice gas, although in these applications the vertex weights are usually not 
chosen as above. 

Another paradigm is coming from the analogy to the self-avoiding walk, 
as in the Brauer model a walk is not permitted to traverse a lattice edge more 
than once. A walk subject to this condition and with no further restrictions 
on the vertices is called a self-avoiding trail (SAT) ^Bj. We will therefore 
refer to the paths in the Brauer model as trails. Here we study the model at 
g = 1 on a bounded domain, and we will be interested in the distributions 
of special points on the boundary that are visited by one of the trails. This 
will be explained in more detail below. 

First we introduce some notation. For given angles a and f3 in the range 
(0, vr) such that a + f3 < n, we define Tq,,/? as the triangle in the upper half of 
the complex plane with vertices at and 1, such that the interior angle at 
is a and the interior angle at 1 is /3. By Wa,i3 we will denote the third vertex 
of Ta^fs- For a given angle cj) £ (0, 7r/2) and integer system size > 0, let 
V = V(f) he the set of vertices {2j cos (p + k exp(i(/)) :j,k^¥\,j + k< N}. 
Properly rescaled, this collection of vertices provides a nice covering of the 
isosceles triangle T^^cf, with base angle cp, see figure ^ 

We define the Brauer model on V with g = 1 as follows. Each vertex 
of V can carry either of the three following configurations of trail segments: 
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Figure 1: Boundary conditions for the self-avoiding trail in a triangle. 

Here, the third vertex state should be interpreted as a crossing. The state 
of each vertex is chosen from {a, b, c} independently from the states of the 
other vertices, with the probabilities for the three states given by 

p(a) = A7r20 p{b) = Xtt {tt - 24>) p{c) = X (p{7r - 24)) (1) 

where A := [n'^ + (p{7r — 2(p)]~^ provides the normalisation. Two distant edges 
may or may not be connected to each other by a trail. The correlations of 
these events are expected to be isotropic in space for large distances, if the 
weights are chosen as in equation and vertices are arranged in space as 
in figure n The anisotropy of the weights is thus precisely compatible with 
the spatial anisotropy of the lattice, see [12j . 

Each configuration of vertex states defines a collection of self-avoiding 
trails on the vertex set V. We will be interested only in those configurations 
in which one of these trails crosses the triangular domain from a given corner 
to the opposite side. To make this trail stay inside the triangle we have to 
impose suitable boundary conditions. We choose the boundary conditions 
as shown in figure ^ Here it is assumed that the system size N is an 
even number, so that each side of the triangle carries an odd number of 
vertices. With these boundary conditions, one trail of the model must enter 
the triangle at the top or the lower-left corner, and stay inside the triangle 
until it leaves at the other of these two corners. 

In our simulations we generate this particular self-avoiding trail dynam- 
ically as follows. Initially, only the vertex states on the boundary are fixed 
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according to the boundary conditions. The vertex states in the interior are 
still undecided. The trail starts from either the top or the lower-left corner 
of the triangle. The steps of the trail follow the trail segments of the vertex 
states. Each time the trail meets a vertex whose state is still undetermined, 
we choose its state according to the probabilities p{a), p{b) and p{c) given 
above, and continue the random walk. The state of this vertex is fixed for- 
ever onwards. We stop the simulation as soon as the trail hits a vertex on 
the side of the triangle opposite the starting point. 

We want to study the distribution of the point where the trail hits the 
side opposite the starting point in the limit when becomes large. In 
fact, we are interested in the distribution of the so-called hull generated 
by the trail up to the stopping time, that is, the collection of points in 
the triangle that are disconnected by the stopped trail from the side of the 
triangle opposite the starting point. This will be motivated in more detail 
in section |2j As we shall explain, the main motivation for our study of 
the hull distribution is the close connection between the Brauer model and 
critical bond percolation, for which the hull distribution is known (although 
rigorously only for critical site percolation on the triangular lattice). 

To make the connection with bond percolation, let us modify the model 
above by setting p{c) = 0, and let p{a) and p{b) be proportional to sin(2(^/3) 
and sin((7r — 2(/))/3), respectively (see for instance ^2 section 5]). For any 
given configuration of the model, one can draw at every vertex either the 
horizontal or the vertical edge connecting the centres of two adjacent rhombi 
(see figure P), such that the drawn edge does not intersect a trail segment at 
that vertex. It is easy to see that the drawn edges constitute a configuration 
of bond percolation on a rectangular lattice together with the dual configu- 
ration. On the boundary we have a percolation cluster along the base and 
right side of the triangle and a dual cluster along the left side. Note that 
as before, there is one special path in the model that crosses the triangle 
between the top and the lower-left corner, and that this path describes the 
interface between the cluster and the dual cluster attached to the boundary. 
This path is called the percolation exploration process. 

Thus we can interpret the Brauer model at g = 1 as a variant of a loop 
model for percolation with the added possibility that the loops may cross. 
Note that these crossings happen only with small probability, since p{c) is the 
smallest of the three weights. One can therefore expect a self-avoiding trail 
to explore space much like the exploration process of critical percolation, 
except that occasionally the trail may cross and possibly re-enter a part of 
space it has already explored. 
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2 Conformal invariance and locality 



In section n we introduced the Brauer model. It is believed that the scaling 
limit of this model is conformally invariant. Moreover, a self-avoiding trail 
of the model has the locality property. The purpose of this section is to 
explain what we precisely mean by these two properties, and to discuss an 
important implication. 

To explain conformal invariance, suppose that for every domain in C 
with continuous boundary and suitable boundary conditions, we are given a 
probability measure on the collection of paths in that domain. Then we say 
that this family of measures (or the family of random paths they describe) 
is conformally invariant if for any conformal map g from a domain D onto a 
domain D' that maps corresponding boundary conditions onto each other, 
the image of the measure on paths in D by 5 coincides with the measure on 
paths in D'. 

Now let T be a triangle in the upper half of the complex plane with 
vertices at 0, 1 and w, and suppose that D is of the form T\A, where A C T 
is such that T \ ^4 is simply connected and w £ {dT \ dA) , see figure Let 
the map gA : T — > T \ ^ fix w and send dT \ [0, 1] onto the boundary arc Ba 
of dT \ dA containing w. Consider a random path W in this triangle which 
starts in w and is stopped as soon as it hits the interval [0,1]. Then we 
say that W has the locality property if the path in T started from w and 
stopped at the first time when it hits dD \ Ba, has the same distribution 
as the path in D started from w and stopped on dD \ Ba- Note that if W 
is conformally invariant, then the latter random path is the image of the 
former by gA- 

For a local, conformally invariant path W in T, let r denote the first 
time when W hits the interval [0,1], and set X := W{t). Then we define 
the hull K generated by the path as the set of points in T that are either 
on Vl^[0, r] or are disconnected from {0,1} by VF[0, r]. Observe that the 
distribution of this hull K is determined if we know for all sets A as in the 
previous paragraph the probability that K <^ A = $. But by the locality 
property, this probability is exactly the probability that the path W va. D = 
T\A stopped when it hits dD \ Ba-, is stopped in the interval [0, 1]. Using 
conformal invariance, this probability equals the probability that X is in 
the interval S';^^([0, 1] \ dA), see figure |2j It follows that for local paths, 
the distribution of the hull K is determined by the distribution of the exit 
point X (see |24) for an illuminating discussion). 

To examine some concrete examples, let T be an equilateral triangle (in 
other words, take w = exp(i7r/3)). Then there are three known local and 
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Figure 2: A hull generated by a local path in a triangle T. The distribution 
of the hull is determined by the probabilities that it avoids sets A such as 
shown on the right. 

conformally invariant paths such that X has the uniform distribution. The 
first example is the trace of Schramm-Lowner Evolution (SLE) for k = 6 
|15| I25j . The second example is a Brownian motion in T which is reflected 
on the left side in the direction given by the vector exp(— ivr/S) and on 
the right side in the direction given by exp(— 2i7r/3) ^^EHj (for details on 
reflected Brownian motion, see section The third example is the scaling 
limit of the exploration process for critical site percolation on the triangular 
lattice in T [25 . Since the distribution of X is the same for these paths, by 
the result of the previous paragraph they generate hulls K with the same 
distribution. 

As we explained in section ^ the Brauer model studied here is closely 
related to critical percolation. One may therefore expect the distribution of 
the hull to be the same for the two models. This can be motivated further 
by the fact that any two local and conformally invariant processes in the 
plane which are started in the origin and stopped upon hitting the unit 
circle generate the same hull, namely that of a stopped Brownian motion 
(see [21] )• In a bounded domain, however, our Monte Carlo study shows that 
the hull distributions for critical percolation and for the Brauer model are 
not the same. Instead, one of the results of this paper is that our numerical 
data are consistent with the hypothesis that in the scaling limit, the hull of 
the Brauer model has the same distribution as that of a Brownian motion 
which is reflected perpendicularly on the boundary of the domain. 

To study the distribution of the hull for the Brauer model, by locality it 
is in principle sufficient to look only at the distribution of the exit point X 
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Figure 3: Definition of the hull K oi a, local path in the triangle T, and of 
the random variables X, Y and Z. 

in a given triangle. However, for obvious reasons we prefer to consider more 
characteristics of the hull. For this purpose we introduce two new random 
variables Y and Z associated with the hull K in the triangle T, as follows. 
We shall denote by \w\Y the distance of the lowest point of the hull on the 
left side of T to the top w, and by \w — 1\Z the distance of the lowest point 
of the hull on the right side to the top w. Thus, all three random variables 
X, Y and Z take values in the range [0, 1]. See figure El for an illustration. 

In section |21 we shall compute the (joint) distributions of X, Y and Z 
for the case of reflected Brownian motion. These distributions can then be 
compared with the corresponding distributions in the Brauer model. The 
results of our Monte Carlo study of the hull distribution are discussed in sec- 
tion^ In section[21we present simulation results for a different distribution, 
the last-visit distribution. 

3 Reflected Brownian motion 

In this section we deflne reflected Brownian motion in triangles and sum- 
marize the main results on these processes. We recall that for given angles 
A and fi in the range (0, vr) such that A + ^ < vr, we define T = Tx^^ as 
the triangle in the upper half of the complex plane with vertices at and 1, 
such that the interior angle at is A and the interior angle at 1 is By 
= w\,fj. we denote the third vertex of 7a,^. 
Let us now introduce the (refiection) angles a,/3 S (0,7r) and set vl '■= 
exp(i(A + a — vr)), vr := exp(— i(/u + /?)). We want to consider a stochastic 
process Z in T which is Brownian motion in the interior, and which is 
refiected instantaneously in the direction given by vl or vr when it hits 
the left or right side of T, respectively. Such a process is called a refiected 
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Brownian motion with reflection vector fields vl and u/j, on the sides of T. 
We write RBM,^^^ to denote this process. Note that the angles a and (3 are 
the reflection angles measured with respect to the boundary. 

More explicitly, an RBM^^^ in T may be defined as follows. Let B be 
standard two-dimensional Brownian motion. Then an RBM^,^^ in T is the 
unique process Z such that 

Z{t) = B{t) + VLYL{t)+VRYR{t), (2) 

where Yl , Yr are continuous increasing real- valued processes adapted to B 
which satisfy Yl{{)) = yR(0) = and increase only when Z is on the left, 
respectively right, side of T |^ equation (2.4)]. This process is well-defined 
up to the first time when Z hits the interval [0,1]. For a characterisa- 
tion and properties of these reflected Brownian motions, see Varadhan and 
Williams (22j. It is straightforward to show conformal invariance and locality 
for these processes, see for instance (23 chapter 5]. 

It was shown by Lawler et al. JH] that if we take an RBM^/3 7^/3 in the 
equilateral triangle 7^/3,77/3, then the exit distribution is uniform. Their 
arguments were generalised to isosceles triangles (i.e. the triangles T\^x) by 
Dubedat [01, and then to the case of a generic triangle (any T\^^ such that 
A + /i < vr, with a natural extension to the case A + /x > vr) by one of 
us These results show that an RBMq,^^ in the triangle Ta^p has the 

uniform exit distribution. By conformal invariance, this also determines the 
exit distribution for an RBMq,^^ in any triangle Ta,^, and using locality we 
can in fact compute the joint distribution of the exit point X and the lowest 
points Y and Z of the hull on the two sides, as we show next. 

It turns out that the distributions of the random variables X, Y and Z 
can all be expressed in terms of conformal transformations of the upper 
half-plane onto triangles T^i^s- By the Schwarz-Christoffel formula PP, the 
unique conformal transformation of the upper half-plane onto T^^s that fixes 
and 1 and maps 00 to w^^s is given by 

F^s{z) = —, ^ . 3) 

By the substitution t 1-^ 1 — u it is easy to show that 

F^^s{z) = l-Fs^^{l-z) and F-}{z) = I - F,^^{1 - z). (4) 

We shall now explain how we can express the distribution functions for X, 
Y and Z in terms of these conformal transformations. 
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Figure 4: This figure illustrates how the joint distribution function of the 
random variables X, Y and Z can be computed. As explained in the text, 
the joint probability P[X < x,Y < y, Z < z] is just g{x) — g{0). 



The idea of the computation of P[X < x,Y < y,Z < z] is illustrated 
in figure 13 Consider an RBM^^^ in the triangle T = Tx,^ started from the 
top w = w\^^. Let a = a{y) and h = h{z) be the points on the left and right 
sides of T at distances \w\y and |u) — l|z from w, respectively. Stop the RBM 
as soon as it hits the counter-clockwise arc from o to 6 on the boundary (the 
thick line in figureE}- Then the probability P[X < x^Y < y, Z < z\ is just 
the probability that this process is stopped in the interval (0, x). 

We now use conformal invariance and locality. Let g = ga{y),b{z) be the 
conformal map of T onto T^^/s that sends a to 0, 6 to 1 and w to 
as illustrated in figure Then we are looking for the probability that an 
RBMq,^^ in Ta^p started from Wa^p and stopped when it hits [0, 1], is stopped 
in the interval {g{0), g{x)). But since the exit distribution of the RBM^^^ is 
uniform in T^^p, this probability is simply g{x) — g{0). Thus, 

P[X <x,Y <y,Z <z]=g{x)-g{0). (5) 

It remains to find an explicit form of the map g = ga{y),b{z)- this point 
it is useful to denote by u the third angle of the triangle Ja,^, that is, 
1/ := vr — A — /u. The explicit form of g can be obtained by suitably combining 
conformal self-maps of the upper half-plane with triangle mappings. How 
this is done exactly is described in figure 121 

By studying this figure we obtain the formula 



P[X <x,Y<y,Z<z] 
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Figure 5: This illustration shows schematically how one obtains an explicit 
form for the map g in terms of the variables y and z. The notations a, b 
and X in the figure are short for F^^{a), F^J^{b) and F^^{x). 



where the images of a and b under F. j can be expressed in terms of y and z 



as 

Sending two of the variables x, y,z to 1 and using the symmetry property Q) 
of the triangle mappings, we obtain 

P[X<x]=F^,^{F^l{x)); (8) 
F[Y<y]=Fp,4F-l{y)); (9) 
F[Z<z]=F^,t,{F-j{z)). (10) 

Observe that these three distributions have particularly simple forms with 
a nice geometric interpretation. For instance, P[Y < y] is just the image 
of y under the transformation that maps the triangle Tj, a onto T/3,0, fixes 
and 1 and takes w,y^x onto Wjs^a- 



4 Hull distribution of the Brauer model 

We now return to the Brauer model introduced in section ^ To compare the 
hull of the model with the hull of reflected Brownian motion, we measure 
the (joint) distributions of the random variables X, Y and Z (introduced in 
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Figure 6: How the random variables X, Y and Z and reflection angles 
a and P are associated to the trails starting from different corners. For 
simplicity we have omitted the normalisation of X, Y and Z. 

section ^ in om' simulations of the Brauer model. The data are collected 
for 10^ independent trails starting from the top and 10^ independent trails 
starting from the lower-left corner on the set of vertices V = for the 
12 different system sizes N = 100, 200, 300, 400, 600, 800, 1200, 1600, 2000, 
2400, 3200, 4000 and 8 different base angles ranging from 10° up to 80° 
with 10° intervals. Below we shall describe more precisely how the data for 
the (joint) distributions of X, Y and Z are collected. 

We distribute the vertices of V on the three sides of the triangle over a 
total of 100 bins, each containing A^/100 vertices. Since there are actually 
+ 1 vertices on each side, this means that one vertex on each side is 
omitted. This is the vertex in the corner of the side where the associated 
random variable X, y or Z is 0. In the simulations we record the number 
of trails Xi that land in the zth bin for each i = 1,2,... ,100. We also 
record the numbers of walks 1^, respectively Zj, such that the ith bin on 
the right, respectively left, side of the triangle (as seen from the starting 
point) is the bin furthest from the starting point which was visited by the 
trail. Figure illustrates how the variables X, Y and Z are associated with 
the hulls of trails starting from different corners. The joint distributions of 
these variables are recorded similarly, but instead of using 100 bins on the 
sides we use 50 bins to reduce memory requirements. 

In this way the simulations build up histograms of the marginal and joint 
distributions of X, Y and Z for the self-avoiding trail that can be compared 
to the corresponding distributions for the reflected Brownian motions. To 
improve the statistics, we first merge together bins of collected data in order 
to distribute the numbers of trails in different bins more evenly. We shall 
explain below how this merging procedure works for the joint distribution 
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Figure 7: The distribution functions P[X < x] (left) and P[y < y] (right) 
for the self-avoiding trail on different triangles with system size = 4000, 
started from the lower-left corner. The solid lines are the corresponding 
distributions for a reflected Brownian motion for which we obtain an optimal 
least-squares fit. See tabled for the corresponding values of a and (3. 

of X, Y and Z. A similar procedure is applied for the joint distributions of 
two of the three random variables and for the marginal distributions. 

For the joint distribution oi X,Y and Z, from our simulations we have a 
total of 50 X 50 X 50 cubical bins that span the unit cube with the variables 
X, Y and Z along the axes. We consider merging together either 10 x 10 x 10, 
or 5x5x5, or 2x2x2 of these cubical bins into larger cubical bins, that 
together with the unmerged bins form a partition of the unit cube. Our aim 
is that each bin in the final partition of the unit cube represents at least 
100 generated trails, 0.01% of the total. Moreover, we want each bin that 
represents at least 1% (10^) of the total number of generated trails to be 
present in the final partition. 

The merging procedure therefore works as follows. First, we consider for 
each of the 10 x 10 x 10 cubes whether it should form a large bin, or be 
built up from smaller bins. So we do the following test: if we would built 
it up from 2x2x2 bins, then would each constituting bin contain at least 
100 trails, or would one of these bins contain at least 10^ trails? If so, the 
10 X 10 X 10 cube is split into 2x2x2 cubes. Then we test whether each of 
these 2x2x2 cubes should form a bin, or be built up from 1x1x1 bins, 
by the same criterion as above. Otherwise, we test whether the 10 x 10 x 10 
cube can be built up from 5x5x5 (rather than 2x2x2) cubes. If so, we test 
whether each of these 5x5x5 cubes should form a bin, or be built up from 
1x1x1 bins, as before. In the end, we achieve a more even distribution of 
trails over a smaller number of bins, as desired. 

As we explained, we want to compare the distribution of the hull gener- 
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a/iT (3/tt 


X 


0.456(1) 0.484(2) 


0.461(2) 0.474(2) 


0.464(2) 0.464(2) 


0.464(2) 0.453(2) 


Y 


0.456(1) 0.446(1) 


0.462(2) 0.445(1) 


0.464(1) 0.438(2) 


0.461(2) 0.431(2) 


Z 


0.447(2) 0.483(2) 


0.445(1) 0.475(2) 


0.439(2) 0.464(2) 


0.431(1) 0.452(1) 



Table 1: Fitted values for a and (3 on different triangles compared for the 
different marginal distributions of X, Y and Z (system size = 4000, trails 
start from the lower-left corner). 

ated by the trails with the hull distribution of reflected Brownian motion, 
in particular in the limit when the system size becomes large. Our working 
hypothesis is therefore that the binned data collected in the simulations is 
predicted by the joint distribution function © for the RBM with parame- 
ters a, f3 for some a and (3 as N ^ oo. To test this hypothesis, we make 
a least-squares fit of the parameters a and (3 for each system size A'^. For 
every triangle T^^,^ on which we have simulated the model and for each of 
the two possible starting points, this gives us a list of the best-fit values for 
a and (3 for 12 different system sizes. We want to investigate whether these 
values of a and (3 converge to a common limit as the system size becomes 
larger and larger. 

The least-squares fits show that the marginal distributions oi X,Y and Z 
for the self-avoiding trail are well described by those of an RBM (see fig- 
ure I7|). More precisely, as the system size gets larger, the value of (the 
difference between the actual and the predicted number of walks in each bin 
squared, divided by the predicted variance in this number, summed over all 
the bins) goes down to a value close to the number of bins, see for instance 
table 121 However, the fitted values we obtain for a and (3 are not constant 
with the system size. More work is therefore needed to test convergence of 
a and /3 to a limit as — >■ oo. 

We also observe from the simulations that the fitted values of a and (3 
for the three distributions of X, Y and Z do not fully agree. It seems that 
the fitted values of a do agree for the distributions of X and Y , whereas 
the fitted values of (3 do agree for the distributions of X and Z. The fitted 
value of a for the distribution of Z agrees with the fitted value of (3 for 
the distribution of Y . See table ^ This observation holds both for trails 
starting from the top of the triangle and for trails starting from the lower- 
left corner at different system sizes. This is another indication that more 
work is required to test convergence of a and (3 as N ^ oo. 

Interestingly, for trails that start in the lower-left corner of an equilateral 
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triangle the fitted values of a and /3 do agree for the exit distribution at 
different system sizes (see table In other words, the exit distribution of 
the trails is symmetric in an equilateral triangle. Note that this result is not 
trivial because our boundary conditions are not symmetric between the left 
side and the base of the triangle, see figure ^ We infer that our choice of 
boundary conditions does not destroy the isotropy of the model, and hence 
does not interfere with conformal invariance in the scaling limit. 

Before we consider the convergence to a limit in more detail, let us also 
consider the joint distributions of X, Y and Z. Since the fits of the marginal 
distributions give different values of a and /? at finite system sizes, it is to 
be expected that the joint distributions measured in the simulations will not 
be well described by those of an RBM. This is indeed what we see when we 
try to fit the data for the joint distributions to the distribution function for 
the RBMs, see for instance table |^ We note that the value of foi' t^is 
fits of the joint distributions is several times the number of bins at small 
system sizes. However, does go down as the system size increases, which 
is a sign that the fits become better for larger systems. 

If the hull of the self-avoiding trail does converge to that of a reflected 
Brownian motion in the scaling limit, then the fitted values of a and /3 should 
converge to a limit value as N ^ oo. From the change in the fitted values 
of a and /? with system size observed in the simulations we can see that 
if there is convergence, then it is very slow. We make the educated guess 
that a and /? converge with the system size as 1/ logj^g ^ (corrections that 
behave as a power of logj^g ^ rather than as a power of itself are consistent 
with earlier findings J7] and with the presence of a zero conformal weight). 

To test convergence, we therefore make a linear fit of the fitted values of a 
and P for the different distributions on different triangles against l/log^g 
We accept the linear fit if it passes each of the following three tests: 

1. We look at the value of for the fit. If this value exceeds the 10% 
probability threshold for the distribution, we reject the fit. 

2. We compare the value of the intercept at = oo from the linear fit 
with the value we obtain if we fit the data to a parabola in 1/ log^Q N. 
If the values do not agree within 1.96 standard deviation (the Gaussian 
5% level) of the linear fit, the fit is rejected. 

3. We do a run test. If the data is well described by a line, then each 
data point should be independently above or below this line with equal 
probabilities. Thus we can predict how many consecutive runs of 
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Table 2: Fitted values for a and /? together with the value of and the 
number of bins for two different distributions of a trail on an equilateral 
triangle, started from the lower-left corner. 



points above the line and below the line to expect. If the probability 
for the number of runs we find is less than 12%, the fit is rejected. 

For more background on these (and other) kinds of tests for data fitting, the 
reader can consult for instance Barlow [2j. 

We have a total of 224 sets of data points for which we attempt a linear 
fit (3 marginal distributions plus 4 joint distributions for the variables X, 
Y and Z, times 8 different base angles, times 2 because the trails can start 
either from the top or the lower-left corner, times 2 variables a. and /?). Of 
these, 39 (17%) give an accepted linear fit in the 12 system sizes, and an 
additional 27 (for a total of 66, i.e. 29%) give an accepted linear fit if we 
leave out the smallest system size. For the other sets of data we can not 
observe the convergence from our simulations without further knowledge of 
how a and /3 should behave as functions of N . 

From each accepted linear fit of q or /? against 1/logj^o-^ obtain a 
value for the intercept with the a- or /3 N ^ oo. Considering only 

the accepted linear fits with all 12 system sizes taken into account, these 
intercept values give a/vr 0.4964(48) and P/ir 0.4952(75) as TV ^ oo. 
If we include the additional accepted linear fits where the smallest system 
size is left out, we obtain a/n 0.4955(54) and P/tt 0.4952(66). These 
results are consistent with the hypothesis that the hull of the self-avoiding 
trail in the scaling limit is the same as that of an RBM^/2,7r/2) that is, a 
reflected Brownian motion with perpendicular reflection on the boundary. 
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Figure 8: Graphs (a), (b) and (c) show the fitted values of /J/vr against 
l/log]^Q for the marginal distributions of X and Z and the joint distri- 
bution of X, Y and Z, respectively, for trails started from the left corner. 
Graph (d) shows the fitted values of q/tt against l/log^o-^ for the joint 
distribution of X, Y and Z for trails started from the top. In those cases 
where a linear fit is accepted, the fitted line is shown as well. 

5 Percolation and the last-visit distribution 

In this section we look at a different distribution for the self-avoiding trail of 
the Brauer model. The motivation for this comes from a further connection 
between the percolation exploration process and ,^^3 discovered by 

Dubedat the analogue of which we want to investigate for the Brauer 
model. But first let us consider the case of percolation. In section ^ we 
explained that the Brauer model studied in this paper becomes a model for 
critical bond percolation if we take the probabilities of the vertex states a 
and h on page 121 proportional to sin(2(/)/3) and sin((7r — 20)/3), respectively, 
and p{c) = 0. The paths defined by the model then correspond to the bound- 
aries of the percolation clusters, or equivalently, to a percolation exploration 
process. The hull generated by a path of the model should therefore have 
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Table 3: Fitted values for a and /? on different triangles for the joint and 
last-visit distributions of the exploration process of critical bond percolation 
at system size = 2000. 

the same distribution as the hull of an RBM.n-/3 j^/s in the scaling limit, as 
we discussed in section |21 

Similar simulations as for the self-avoiding trails allow us to test this 
hypothesis. We have measured in the same way as before the joint distri- 
bution of the variables X, Y and Z on different triangles with fixed system 
size N = 2000. Our hypothesis is that this joint distribution is predicted by 
equation © for reflected Brownian motion with reflection angles a and /3, 
where a and /3 should be equal to vr/S. Table 01 shows the results of a least- 
squares fit for the reflection angles a and /3, which agree with the hypothesis 
that the hull distribution is the same as that of an ,^^3 . 

It was shown by Dubedat jjj that there exists a rather surprising further 
connection between reflected Brownian motion and the exploration process 
of critical percolation. To explain this connection, suppose that is a local 
process in the triangle T = started from w = wx^^, where A,// are fixed 
angles such that A + /i < vr. Let r := mf{t > : W{t) G [0, 1]}, and let a be 
the last time before r when W visits the boundary of T. Let E denote the 
event that W{a) is on the right side of the triangle. Then we can consider 
the probability of the event E conditioned on {VF(r) = x}. We call this 
conditional probability, which is a function of x, the last-visit distribution 
of W in T. 

For reflected Brownian motion this last-visit distribution can be com- 
puted using the fact that the exit distribution is uniform in a well-chosen 
triangle. For the case of an RBM^/3^/3 in an equilateral triangle, this com- 
putation was done by Dubedat in The generalization to an RBMq^^ in 
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Figure 9: Fitted values of a/vr (left) and P/tt (right) plotted against 
1/ log]^Q N for the last-visit distribution of trails starting from the lower-left 
corner of different triangles. 



the triangle Tx^^ started from wx^n is straightforward, and gives 

P[E I Wir) = x] = F^_a,._;3 {P^lix)) , (11) 

where F^^g, given by equation (jH)), is the unique conformal transformation 
of the upper half-plane onto ,5 which fixes and 1 and sends cxd to w^^s- 

The last-visit distribution can also be computed rigorously for chordal 
SLE(6), as was shown by Dubedat f^. The obtained formula is exactly the 
formula (|11|1 in the case of an RBM^^^ .^/^ (i.e. a = (3 = ir/S). This fact is 
rather surprising if one considers the very different ways in which the two 
processes explore space: whereas an SLE(6) process never crosses itself, an 
,r/3 crosses itself many times. In particular, it is known that on 
the event E, the last point visited on the right side of the triangle T by 
an SLE(6) process must be also the lowest point of the hull on the right 
side. For an however, the last point visited on the right side is 

almost surely not the lowest point of the hull. 

Let us now return to our simulations of the exploration path of critical 
percolation. It is quite easy to count in these simulations the number of 
paths that land in a given bin and visited the side to the left (and not 
to the right) of the starting point last before reaching the opposite side. 
Divided by the total number of paths that land in a given bin, this binned 
data can be compared to the last- visit distribution of equation (|11|) . Note, 
however, that this binned data represents the conditional probability of E 
given that the path lands in an interval rather than at a given point x. To 
make the comparison we therefore have to integrate (|11() with respect to the 
exit distribution of the exploration path. 
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We take the exit distribution for the paths in the simulations to be well 
described by the least-squares fit for the marginal distribution of X. Then we 
can integrate numerically (with a and (3 as parameters) with respect 
to this fixed exit distribution, and obtain from this the probability of E 
conditioned on the event that the path lands in a given bin. This can be 
used to make a least-squares fit for the parameters a and [3 of the last-visit 
distribution. Table 01 shows that our simulations for critical percolation are 
in agreement with the fact that the last-visit distribution is given by that 
of an 

A similar analysis can be carried out for the self-avoiding trails of the 
Brauer model. As before, here we face the problem that the model converges 
only slowly. We therefore make least-squares fits of a and jS for different 
system sizes, and plot the results against l/log^oA^, see figure IHl It is 
clear that here we can not conclude what the scaling limit is, without an 
analytic prediction of how the parameters a and (5 behave as functions 
of N . However, our results do not exclude the possibility that the last-visit 
distribution for the self-avoiding trails in the scaling limit is the same as 
that of an RBM^/2,7r/2) just like the distribution of the hull. 

Conclusions We have carried out a classical Monte Carlo study of the hull 
distribution for the Brauer model at g = 1 . Our results are in agreement with 
the hypothesis that the hull distribution coincides with that of a refiected 
Brownian motion with perpendicular refiection. We have tested our methods 
for percolation, and found good agreement with the fact that in this case the 
hull must have the same distribution as that of a Brownian motion which 
is reflected at angles of 60° with respect to the boundary. We have also 
studied the last-visit distribution, but for the Brauer model our simulation 
results are inconclusive. 
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